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■ Computation of collisional energy loss in a finite size QCD medium has become crucial to obtain 
I reliable predictions for jet quenching in ultra-relativistic heavy ion collisions. We here compute this 

energy loss up to the zeroth order in opacity. Our approach consistently treats both soft and hard 

I contributions to the collisional energy loss. Consequently, it removes the unphysical energy gain 

■ in a region of lower momenta obtained by previous computations. Most importantly, we show 
I that for characteristic QCD medium scales, finite size effects on the collisional energy loss are not 
• significant. 

m ■ 

o : 

^ ; 1 Introduction 

1^ ' The suppression pattern of high transverse momentum hadrons is a powerful tool to map out the 

' ^ ■ density of a QCD plasma created in ultra-relativistic heavy ion collisions (URHIC) ^-0- This sup- 

pression (called jet quenching) is considered to be mainly a consequence of medium induced radiative 
^ I energy loss of high energy partons |1|-[Z|. However, recent non-photonic single electron data jHl El 

^ ' (which present an indirect probe of heavy quark energy loss) showed large disagreement with the ra- 

• i-H . diative energy loss predictions, as long as realistic values of parameters are assumed |1U| . This raised 

^ I the question of what is the cause for the observed discrepancy. 

d • Recent studies ^] suggested that one of the basic assumptions that pQCD collisional energy 

loss is negligible compared to radiative jl3j may be incorrect. In |lllll2j it was shown that, for a range 
of parameters relevant for RHIC, radiative and collisional energy losses for heavy quarks were in fact 
comparable to each other, and therefore collisional energy loss can not be neglected in the computation 
of jet quenching. This result came as a surprise because from the earlier estimates |18j . the typical 
collisional energy loss was erroneously considered to be small compared to the radiative energy loss. 
In it was consequently suggested that the inclusion of collisional energy loss may help in solving 
the "single electron puzzle". However, in that study (as well as jllj-|18j) finite size effects were not 
taken into account. 

A recent paper by Peigne et al. j^Oj is the first study that made an attempt to include finite size 
effects in the collisional energy loss. This work suggested that collisional energy loss is large only in 
the ideal infinite medium case, while finite size effects lead to a significant reduction of the collisional 
energy loss. However, this paper did not completely separate collisional and radiative energy loss 
effects. Consequently, it remained unclear how important are the finite size effects on the collisional 
energy loss. 

Therefore, it became necessary to consistently compute (only) the collisional energy loss in finite 
size QCD medium. Additionally, this paper aims to address whether -and to what extent- there is an 
over-counting between collisional and radiative energy loss computations. 
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The outline of the paper is as fohows: In Section 2, we will compute the collisional energy loss in a 
finite size QCD medium. In Section 3, we will consider the special case when a particle is produced at 
xq = — oo (infinite medium case). We will show that in special limits, our calculations recover previous 
results |141ll5j . However, contrary to jl4l I15j . our computation does not encounter unphysical energy 
gain in the low momentum region jl4l I15j . In Section 4 we will give a numerical study of the collisional 
energy loss in both finite and infinite QCD medium. Contrary to the results obtained in Peigne et 
al. |2nj we will show that finite size effects do not have a significant effect on the collisional energy 
loss. The conclusions and outlook are given in Section 5. 

2 Collisional energy loss in finite size QCD medium 

In this Section we will compute the collisional energy loss (up to the 0*'^ order in opacity) when the 
jet is produced in a finite size dielectric medium. The contribution to this energy loss comes from one 
Hard-Thermal Loop (HTL) gluon propagator (see Appendix A), which is the reason why we call it the 
Qi/i order in opacity collisional energy loss (note the analogy with the 0*^ order in opacity radiative 
energy loss |2^-(2S], which is further discussed in Appendix A). 

In this computation we use the most intuitive approach, i.e. we compute the diagram |Me;| shown 
in Fig. Note that the "blob" represents the effective gluon propagator. A proof of the validity of 
this approach is given in Appendix A. This approach has already been used without proof in 1^ . 
under the justification that it reproduces the same results as the imaginary time formalism. 
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Figure 1: Feynman diagram for the amplitude that contributes to the collisional energy loss in QCD 
medium. The large dashed circle ("blob") represents the effective gluon propagator j21j . 

Similarly as in [SHI we introduce the finite size medium by starting from the approach described 
in We consider a static medium of size L, and assume that collisional energy loss can occur 

only inside the medium. The Feynman diagram \Mei\ (see Fig. ^ then represents the source J, which 
at time xq produces an off-shell jet with momentum p, and subsequently (at xi > xq) exchanges a 
virtual gluon with parton in the medium with momentum k. The jet and the medium parton emerge 
with momentum p' and k' respectively. Since our focus is on heavy quark jets with mass M, we here 
neglect the thermal shifts of the heavy quark mass. 

We assume, as in j^, that J changes slowly with p', i.e. that J(j)' + g) ~ J{p)- Since we consider 
both soft and hard contributions, we take into account spin effects. The computation that we present 
in this paper is gauge invariant |21], but for simplicity we further use Coulomb gauge. 

The effective gluon propagator shown in Fig. 1 has both transverse and longitudinal contribu- 
tions |28j-j32j. In Coulomb gauge the gluon propagator has the fallowing form: 

D^^'-iuA) = -P^'^Ar(^,q) -Q^'^AL(u;,q), (1) 
where q = (to, q) is the 4-momentum of the gluon, while At and Al are effective transverse and 
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longitudinal gluon propagators given by P^ : 



^ ^2 2q2 ^ 2|q| 'u; + |q|'^' ^ ^ 



{u}"^ — a^)y? ^ uo ^^lO;— |q| 

|q| 
|q| 



A^^=q2 + ^2(l + -r-ln|^— ^1), (3) 
2|q| uj + |q| 

where = g'^T'^{l + ^) is the Debye mass. 

The only nonzero terms in transverse {Pf_Ly) and longitudinal (Q^u) projectors are the following: 

= 1. (5) 

The matrix element for this 0*'^ order in opacity collisional process can then be written in the 
following form (for simplicity we here omit color factors, whose contribution we will add in the end) 

iM,i = j d^xo J(xo) d^xi i-i) J 0(ii - to) e{L/v - {h - to)) 

u{k', A')e*^'^2 igY u{k, \)e-'^''^ . (6) 

Here p, s, k and A are the 4-momenta and spins of the incoming jet and medium parton, while the 
corresponding primed variables correspond to outgoing jet and medium parton (the medium parton can 
be quark, antiquark or gluon). The medium partons are considered to be massless, i.e. 4- momentum 
k (k') is assumed to be A: = (|k|,k) (A;' = (|k'|,k')). The condition that the interaction between jet 
and medium parton has to occur inside the QCD medium of finite size L is implemented through 
the second 9 function giving maximal interaction time of (ti — to)max = L/v. We will further define 
X = Xi — xq = {t, x). 

The Eq. © simplifies to 

{2Ti)H{k' -k-q) D^,{q) u{p\ s'^^^u^p, s)u{k' , X'^'u^k, A) 

1 1 _ p-i{E-E'-uj)L/v 

where E = \l M"^ + p^, M is the jet mass, p = p' — (k' — k) and oj = |k'| — |k|. 
In the last step we used J{p' + g) ~ J{p') |2Z| and defined 

M = g^D^^{k' - k)u{p', s')-f''u{p, s)u{k', X')Yu{k, A). (8) 

In this paper we consider the case of highly energetic jets, where |q| <^ E. In this limit E' becomes 
E' ~ E — V ■ q. Here v is the velocity of the initial jet, i.e. the jet 4-momentum p is equal to 

p_ ( M Mv \ 
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Further, the matrix element given in Eq. ^ has to be squared, averaged over initial spin s of the 
jet and summed over all other spins. 



V • qj 

spms spins 



where | Yl spins l-^P is given in Appendix B. 

The differential energy loss is equal to dE^i = tu dT^i , where collisional interaction rate dT^i can be 
extracted from Eq. © as (see [77] ) 

2^' ' 27r 32£^' (27r)^2A; 27r 32A;' ^ ^"^^ ^"^^ ^ ^ 

Here 

<'''v^ = <'«i-Ap')i^p$5i5^. (11) 

and d^ = 3 (for three dimensional representation of the quarks). In Eq. (jlOj) . n\q{k) is the equilibrium 
momentum distribution at temperature T of the incoming parton ^, where ^ denotes quark, antiquark 
or gluon. (1 ± niq{k')) is a factor associated with the outgoing parton, where + corresponds to the 
gluon contribution, and — to (anti)quark contribution. 

In this paper, we are interested only in the computation of the collisional energy loss. It is 
straightforward to show that in the collisional energy loss calculations, the ^n\q{k') part in (l±nfg(fc')) 
can be dropped because the corresponding term in the energy loss integrand is odd under the exchange 
of k and k' , and integrates to zero (see also 24 ) . Therefore, for the purpose of computing the collisional 
energy loss, we can replace (1 it n\q{k')) by 1 in Eq. (|Tn|) . which leads to 

d^NjdVel ~ - > MeiP- r|- rr^ -.^ Ueaik). (12) 

2 ^ ' {2^Y2E' {2^Y2k {2^f2k' ' ^ ^ 

spms ^ ' ^ ' ^ ' 

Here, neq{k) = Y(^=q q gfAq{k) is the equilibrium momentum distribution [21] at temperature T in- 
cluding quark, antiquark and gluon contributions (see Eq. (PUjl l. 



The collisional energy loss can now be obtained from Eqs. Q, pij> and (|12j) . leading to 

^ I f d^k. /■ d^k' sin[(a;- v-q)T^]2 l 

E^ J 27r)32A: 'J (27r^2A;' u; - v • q)2 2 ' ^ ' 



spins 



Note that in Eq. ((T^ we added a color factor Cr and that q = k' — k. 

Equation (|13|1 can be further simplified by noting that, in a static medium, the collisional energy 
loss does not depend on the direction of v. After evaluating the ^ ^ |A^|^ and averaging the integrand 
over the directions of v, we obtain (see Appendix B) 

AE^l = I neg(|k|)(i|k| / |q|(i|q| / u)duj + / |q|d|q| / ^^duj 

Jo \Jo J~\q\ J|k| -^Iql-slkl / 

(l A.(,)|^M±^^^, + |^,(,)|2 (|qP-.-)((2|kK.)2 + |qP) _ ^ ^^^^ _ 



(14) 
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where J2 and are given by Eqs. H34II36() in Appendix B, and |q|maa; is given by the fohowing 
formula 



|4|_ = M.[E,W^|. (15) 

l-v + 2\k\/E 

Further numerical study of the collisional energy loss in a finite size QCD medium is given in 
Section 4. 



3 Collisional energy loss in infinite QCD medium 

Previous calculations of the collisional energy loss, in particular jlll 1141 [T3] were done for an infinite 
QCD medium. A problem with these calculations was that they produce unphysical energy gain 
in the lower momentum regions. Additionally, these computations leaded to different results and 
consequently provided quite a large uncertainty in the heavy quark (especially bottom) collisional 
energy loss. 

In this Section and Appendix B.l we present an improved calculation of collisional energy loss 
in an infinite QCD medium, with the goal of 1) removing the problems associated with previous 
calculations and 2) producing reliable infinite medium results which we will in Section 4 compare with 
the collisional energy loss in a finite medium. 

In the case of an infinite QCD medium, the collisional energy loss per unit length is computed 
by assuming that the jet is produced at xq = —00. The energy loss for a finite size medium is than 
(simplistically) calculated by multiplying this with the thickness L of the medium. 

In Appendix B.l we present an improved calculation of the collisional energy loss per unit length 
in an infinite QCD medium. The following result is obtained: 



^ / neg(|k|)d|k| / d\q\ / cuduj + / d\q\ / coduj 

Jo \Jq J-v\q\ J\ii\/(l+v) J|q|-2|k| 



dL bV^TT-" Jq yjQ J-v\4\ J\Il\/(1+v) ^|q|-2|k| 

/ ,„ (2|k|+^)^-|qp , ,^ ^„ (|qp-a;2)((2|k|+^)2 + |qp) ^ \ 
\^L{q)\ ^ + \^Tiq)\ 7^7 [v |q| - UJ ) ] . (16) 



To compare our result with the computations done in jl4l I15j . we here introduce an arbitrary 
intermediate momentum scale |q|* jT^]. The contribution from |q| < |q|* is denoted soft, while 
contribution from |q| > |q|* is denoted hard |15j . 

The soft contribution is given by 

r°° - - /■I'll* /■''I'll 



'^el ^ _9_ 

dL 6 v"^ vr^ 



/ neg(|k|)ci|k| / d\^\ I ujdio 

Jo Jo J-v\^\ 



„2 (2|k|+a;)^-|qp , ,^ ^„ (|qp - ^^)((2|k| + o;)^ + |qp) ^ \ 
\^Liq)\ + |AT(g)| ^^j4 [v |q| -UJ )\ .(17) 



while the hard contribution is given by 
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^' — ' neg(|k|)d|k| ( / d|q| / udu + I (i|q| / udu 

\J|q|* J-v\ci\ J\^\ ^|q|-2|k| 

The soft contribution can be further simphfied by keeping only the even contributions in the ui 
integral (the odd contributions vanish under symmetric integration) 



dL 3 v'^ vr^ 



k|neg(|k|)d|k| d\q\ 
Jo 





.^dc (|Al(.)P + ^(i-^)(.^-^ 



)|AT(g)| 



n2 /-Iql* /""Iql / 1 , ,2 , ,2 

^ /.I, / d|q| / ..2^^ f |^^(^)|2 + _ ^)(^2 _ ^)|AHg)P ) (19) 



OTTU Jo -^--ulql V ^ 1^1 111 

where in the last step we used the fact that 

9rri9 TV, 
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Equation H19|) agrees with the Eq. (55) in 24_. 



k|ne,(|k|)d|k| = ^^(1 + ^). (20) 



For the hard contribution we use that, in the limit of large momentum transfer (|q| ^ |q|*)) 
|AL(g)| —>■ j^p and |A7-(g)| — > ^^^^j^- It is than straightforward to show that the hard contribution 
reduces to the Eq. (17) in i.e. 

= / "e<? k d k / ^ / uoduj + I ^ uduj \ 

dL Gv^tt'^Jq ^1^1* |qr j-^|q| 72|k|/(i+i,) |qr J|qi-2|k| J 

fi44_i_^ w ,3Ma±^_(i_oMw±^V (21) 

\4|q|2 4 2 |q|2-a;2 |q|2 |q|2-u;2 i 

Equations (17) and (55) from j24j were used in Ref. |15j to obtain their Eqs. (8) and (12). That is, 
while our Eq. (|T6|) is more general, in special cases (i.e. Eqs. (fT9|) and (|2T|) ) it reproduces results 
from [15:- 

The computation in considered only the soft gluon limit, and replaced |q|* by |q|maa;) where 
|q|»nax is given by Eq. ((T3|) . Consequently, for the purpose of comparison with ^1], we replaced |q|* 
by |q|max in Eq. 1)19^ . Additionally, the problem with this approach is that, in the high momentum 
|q| region, the method ^3] is not able to treat the lower oj bound properly (compare Eq. ()19p with 
Eq. (|16() where lo bounds are properly treated). To overcome this problem, the calculation in ^Ij was 
limited to the forward emission only (i.e. uj > 0). If this limit is also taken into account, our Eq. (|19|1 
reproduces Eq.(4.1) in [Tl] . 

In summary, in the known limiting cases, our result (i.e. Eq. ()16() ') reduces to the results published 
in |14| II 5j . The advantage of our result over [Tlj is that it includes the hard contribution and consis- 
tently treats the integration limits. Comparing to JH]) the advantage of our result is that it does not 
make a sharp transition from soft to hard limits, and consequently it does not require the introduction 
of an unphysical momentum scale |q|* as in jl5j . 
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4 Numerical results 



In this section we give a numerical study of the colhsional energy loss in a QCD medium as presented 
in Sections 3 and 4. To do this, we further assume that the QCD plasma is characterized by T = 
0.225 GeV, Nf = 2.5 and a = 0.3. For the light quark mass we take M = ^£)/\/6, where fJ-n = 
g^T^(l + ^) « 0.5 GeV is the Debye mass. The charm mass is taken to be M = 1.2 GeV, and the 
bottom mass is M = 4.75 GeV. 



4.1 Collisional energy loss in infinite QCD medium 
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Figure 2: Fractional collisional energy loss is shown as a function of momentum for light, charm and 
bottom quark jets (left, center and right panels respectively). Dash-dotted curves are obtained by 
using Eq. (|16|) from this paper. Dashed curves correspond to Eqs. (8) and (12) from JH]) while dotted 
curves are obtained by using Eq. (4.1) from jjlj. Assumed thickness of the medium is L = 5 fm. 



In Fig. [21 we compare our collisional energy loss results in an infinite QCD medium (Eq. (|16|) ) 
with previous computations by ^^E]- We see that, while both BT [TSj and TG computations 
lead to unphysical negative energy loss results in the low momentum region, our computations give 
positive collisional energy loss in the whole jet momentum range. This is particularly important in 
the bottom quark case where the unphysical behavior persists up to 5 GeV in BT case and 
up to 2 GeV in TG jl4j case. The reason for this behavior is that only the leading logarithmic 
contribution was considered in the final steps of both BT and TG computations. Note that the problem 
of unphysical energy gain was addressed in Ref. by including fully dressed gluon propagator in 
their calculations. However, that approach leaded to another problem, since the unphysical momentum 
scale |q|* appeared in the collisional energy loss results |33j . Therefore, our results present a first 
complete solution to the unphysical energy gain problem. 

Our numerical results agree with BT only in the asymptotic regions, which is likely the consequence 
of the following: 1) BT made a sharp (instead of continuous) transition from soft to hard limit and 2) 
they introduced a sharp boundary in the energy loss computations depending on whether the initial 
jet energy is much larger /smaller than M^/T. 

Despite the fact that the BT computations are more up to date and treat the collisional energy loss 
more consistently than TG, we see that our results show better agreement with TG computations. 
Particularly, in case of light and charm quark jets, there is a quite good agreement between our 
results and those of TG Jl]. The good agreement is probably because the forward emission only (see 
Section 3) provides a quite plausible estimate for the collisional energy loss. However, for bottom 
quarks we see that neither BT nor TG computations provide a reasonable estimate for the collisional 
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energy loss. Therefore, in this case, the more accurate computation of cohisional energy loss (i.e. our 
Eq. (dl)) is needed. 

4.2 Collisional energy loss in a finite QCD medium 



0.4 f 




Figure 3: Fractional collisional energy loss is shown as a function of momentum for light, charm 
and bottom quark jets (left, center and right panels respectively). Full curves correspond to finite 
medium case (see Eq. dJ), while dash-dotted curves correspond to infinite medium case (see Eq. (jTBI) ). 
Assumed thickness of the medium is L = 5 fm. 
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Figure 4: Fractional collisional energy loss is shown as a function of thickness of the medium for 
light, charm and bottom quark jets (upper, middle and lower set of curves respectively). Full curves 
correspond to finite medium case (see Eq. (jl4j) ). while dash-dotted curves correspond to infinite 
medium case (see Eq. ()16() ). Initial momentum of the jet is 10 GeV. 

Figures El and |3 show the comparison between collisional energy loss in infinite and finite size 
QCD medium. Contrary to [20], we find that a finite size medium does not have a large effect on the 
collisional energy loss. The discrepancy between our results and those presented in [20] is due to the 
fact that what was called collisional energy loss in [201, ^^^^ combination of collisional and part 
of the 0*^ order radiative energy loss. Actually, the calculation in 20 does not present a complete 0*'' 
order energy loss either, since transition radiation [231 was not included in their computation. 
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Figure 5: Fractional collisional energy loss is shown as a function of thickness of the medium for 
charm quark jet. Dash-dotted curves correspond to infinite medium case (see Eq. (|16|)). while dashed 
curve show what would be the collisional energy loss in finite size medium if term 2ixiJ2 — = ^■ 
Dotted curve shows the extra contribution to the collisional energy loss in finite size medium, due the 
fact that 2bjJ2 — Jzi^^- Initial momentum of the jet is 20 GeV. 



Contrary to naive expectations, from Figs. El and 0] we found that collisional energy loss in a finite 
size medium can be somewhat larger than in an infinite medium. The reason is that in Eq. there 
exists a term 2a; ^2 — Jz- If this term were equal to zero (as in the case of infinite medium), the 
energy loss in a finite medium case would always be smaller than in an infinite medium, as naively 
expected (compare dashed and dot-dashed curves in Fig. ISJ. However, in the finite medium case, the 
term 2ijjJ2 — giving a noticeable positive contribution (see the dotted curve in Fig. El) which 

will lead to somewhat larger (overall) energy loss in the finite medium case. 

To further discuss this, let us look at the Eqs. (|34I.S6j) (Appendix A) in the finite medium case, and 
compare them to the corresponding Eqs. (|52I54|) in the infinite medium. The 5 function in Eqs. ()52I54|) 
ensures energy conservation, which is satisfied when the jet is produced at —00. Consequently, in this 
case 20JJ2 — iTs = (see Eqs. andjM})- However, when the jet is produced at finite time xq, time 
translation invariance is broken, and therefore the energy of the collisional process is not conserved, 
leading to 2ujJ2 — Jzi^^. 



4.3 Comparison between collisional and radiative energy loss in a finite size QCD 
medium 

In Appendix A we showed how to separate the contributions to the collisional and radiative energy 
loss. In this section we directly compare these two contributions in the case of finite size QCD medium. 

To compute the net radiative energy loss, we note that there are three important effects that 
control this energy loss in a QCD medium. These effects are the Ter-Mikayelian effect j21j . transition 
radiation |521 and medium induced radiation |34j . In |531, we combined these effects to compute 
the net radiative energy loss. We here use these results for the purpose of further comparison with 
the collisional energy loss. Note that in these computations, in order to simulate confinement in the 
vacuum, we introduce an effective gluon mass in the vacuum m^^^ « ^qcd = 0.2 GeV (for more 
details see |2I])- 
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Figure 6: The comparison between collisional and radiative fractional energy loss is shown as a func- 
tion of momentum for light, charm and bottom quark jets (left, center and right panels respectively). 
Full curves show the collisional energy loss, while dot-dashed curves show the net radiative energy 
loss. Assumed thickness of the medium is -L = 5 fm and A = 1.2 fm |19j . 
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Figure 7: The comparison between collisional and radiative fractional energy loss is shown as a 
function of the thickness of the medium. Light, charm and bottom quark cases are shown on the left, 
center and right panels respectively. Full curves show the collisional energy loss, while dot-dashed 
curves show the net radiative energy loss. Mean free path is A = 1.2 fm ^U]. Initial momentum of the 
jet is E = 10 GeV. 

In Figures ini and [7| we show the collisional and radiative energy loss as a function of jet energy and 
thickness of the medium, respectively. We see that collisional energy loss is comparable with the net 
radiative energy loss in the medium. Therefore, the collisional energy loss contribution is significant 
and must be included in the computation of jet quenching in a QCD medium. 

In particular, we note that in the lower momentum (i.e. p < 10 GeV) range, the collisional energy 
loss dominates the radiative one. At RHIC, jet suppression is mostly measured in this (lower) mo- 
mentum range. Therefore, contrary to previous studies |3|-[Z|) our results indicate that it is collisional 
instead of radiative energy loss which gives the dominant contribution to the observed suppression 
values. 
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5 Conclusion 



This paper addressed the order contribution to the colhsional energy loss in a finite size QCD 
medium. The interest in the colhsional energy loss has been renewed by the recent studies jllL \V2\ . 
particularly in the context of the single electron puzzle [5]-jlUj. In Refs. I12j it was claimed 
that, contrary to the previous beliefs, for the parameter ranges relevant in URHIC, radiative and 
colhsional energy loss become comparable. However, a recent study by Peigne et al. i20j suggested 
that colhsional energy loss is large in the ideal infinite medium case, while the finite size medium 
effects lead to significant reduction of the colhsional energy loss. The paper [20], however, did not 
completely separate colhsional from radiative energy loss. 

Additionally, even in the infinite medium case, the problem of colhsional energy loss was not 
completely solved. Previous computations obtained unphysical results in the low momentum re- 
gions [HI [inj ; and an approach to solve this problem [3S| leaded to the results dependent on the 
unphysical momentum scales. In addition, these computations introduced quite a large uncertainty in 
the heavy quark (especially bottom) colhsional energy loss, since they leaded to noticeably different 
numerical results. 

The above reasons and the previously discussed single electron puzzle, motivated us to provide 
a detailed study of the 0*'' order colhsional energy loss in a finite size QCD medium created in 
URHIC. First, in Appendix A we showed that, though O*'^ order colhsional and radiative energy loss 
contributions come from the same one-loop HTL diagram, there is no overlap between colhsional and 
radiative energy loss computations. More specifically, while O*'^ order colhsional energy loss comes from 
the processes which have the same number of incoming and outgoing particles, the radiative energy 
loss has one gluon more as an outcome of the process. Additionally, we showed that in the O*'^ order 
calculations, there are no interference effects between colhsional and radiative energy loss, which is 
different from a result in the recent paper [3^. The absence of interference effects comes from the fact 
that, contrary to [HEl; in our study we consistently treat the gluon dispersion relation in the medium. 
This leads to the following conditions: 1) for the O*'^ order colhsional energy loss contributions, the 
energy of the exchanged (virtual) gluon has to be smaller, or equal, to the gluon momentum, and 2) 
for the radiative energy loss contributions the energy of the radiated gluon has to be larger than its 
momentum. Therefore, these two contributions take non-zero values in non-overlapping regions, and 
consequently cannot interfere with each other. 

In the case of infinite medium, our computation interpolates smoothly between soft to hard contri- 
butions and, contrary to does not require the introduction of an arbitrary intermediate momentum 
scale. Additionally, our computation treats the lower momentum region consistently, removing the 
unphysical energy gain results obtained in previous computations |14[ II 5j . 

In the case of finite size QCD medium, contrary to the study by Peigne et al. [20] we showed 
that finite size effects have small effect on the colhsional energy loss. Therefore, consistently with the 
claims in Refs. 112j and our recent single electron suppression estimates |19| . we here showed that 
colhsional energy loss is important, and has to be included in the computation of jet quenching. 
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A HTL contribution to the collisional energy loss 

In this section we will derive the formula for the lowest order collisional interaction rate. The zeroth 
order contribution to both radiative and collisional rates comes from the diagram M given in Fig. |S1 
We will below use this diagram as a starting point to separate contributions of collisional and radiative 
energy loss. 




J(P) p' J'(P) 

Figure 8: One Hard Thermal Loop (HTL) diagram. 



Diagram M corresponds to where M„ is the amplitude of the diagram shown in Fig. |UJ 




Figure 9: Diagram M„ 
The definition of "black circles" in Fig. |21 is shown in Fig. IIUI 




Figure 10: Definition of the "black circles" in Diagram M, 



The diagram M contains both collisional and radiative order contribution to the jet energy 
loss. It is useful to look at the simple n = 1 case (see Fig. [TT|) to better understand this. 
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k^k kk'k kk'k kk 




Figure 11: n = 1 contribution to the HTL diagram. 



The contribution to the colhsional energy loss is obtained by "cutting" (i.e. putting on shell) the 
propagators of parton k' and p' . On the other hand, the radiative contribution is obtained by putting 
the parton propagator p' and the gluon propagator q on shell. From this, it follows that collisional and 
radiative contributions come from different diagrams. Furthermore, from the conservation of energy 
and momentum it can be shown that cutting the gluon propagator q, leads to the condition |u;| > |q|, 
while cutting the propagator of parton k' leads to the condition |ti;| < |q|. Consequently, there is no 
interference (and over-counting) between collisional and radiative contributions^. The computation 
of the radiative 0*'' order energy loss has already been a subject of our previous work [231 121j . So, 
the contributions from the diagrams which give raise to the radiative energy loss, will not be further 
addressed here. 

As we can see from the right side of the Fig. IIH there are two contributions from diagram Mi to 
the collisional rate. These two contributions can be combined into one by allowing that the energy 
of the gluon can take both positive and negative values. Therefore, the contribution to the collisional 
rate from diagram Mi (cITmi) is equal to 

where Meq is the Feynman diagram shown in Fig. 1121 and uj G (— |q|) |q|)- Note that in the above 
equation each of the terms in neq{k) = '^^=q q gniq{k) should be multiplied by an extra factor (1 it 

'rAq{k')) for the outgoing medium parton (see Section 2). Here ^ corresponds to anti(quark) or gluon. 
The + sign is associated with gluon and the — sign with anti(quark) contributions. However, the 
ibn|g(A;') term in (1 ± n\q{k')) does not contribute to the collisional energy loss (for more details see 
Section 2 and Ref. |21])- Therefore, in Eq. (|22|) we keep only the contributions that give raise to the 
collisional energy loss. 




Figure 12: Feynman diagram for the lowest order collisional energy loss. 
In the same way, it can be shown that the contribution to the collisional energy loss from diagram 



^Note that in this paper, we treat only the O"* order contribution to collisional energy loss. It is, however, possible 
that interference effects and/or over-counting between collisional and radiative energy loss contributions would occur in 
the higher order computations. Higher order contributions are a separate non-trivial problem, which is not considered 
in this paper. 
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Mn is equal to 



n-l 



{27r)^2E' (27r)32fc (27r)32/c' 



(23) 



i=0 



here the diagram Me„ is shown in Fig. 1131 

k 



k' 




q /' 




Figure 13: Feynman diagram for the coUisional energy loss with n interactions with the medium. 
Since M = Yln'=o -^^i' contribution to the collisional energy loss from the diagram M is equal 



to 



d^Nj dV 



d^p' d^k d^k' 



oo n— 1 



(2^)32^' (27r)32A; {2TTf2k 
d^p' d^k d^k' 



n=0 i=0 

oo n 



{2tt)^2E' {2'Kf2k {2Tif2k 



We next want to prove that d^Nj dV = ^2ff'2E' (2^)^2fc (2^)^2fc' '^e? 



(24) 



n=0 i=0 

{k) |Me/p, where (see Fig. ITU) 



Mez = Me„ 



(25) 



n=0 




Figure 14: Feynman diagram M^i for the collisional energy loss in QCD medium. The large dashed 
circle ("blob") represents the effective gluon propagator fIT\ . 



To prove the above, we will first compute |Me;|^ 



i^^'i' = E E ^^M, = E E 

i=0 j=0 



(26) 



i=0 n=i 
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where in the last step we defined n = i + j — >j = n — i. 
Since 



oo oo oo n 



EE = EE. 

i=0 n=i n=0 i=0 

we can conclude that 

j3_r/ J3it ^3i7/ °° 
(27r)32£;' (27r)32A; (27r)32A:' ^ ^ 

which is what we wanted to prove. Therefore, the collisional interaction rate can be obtained by 
an intuitive approach of computing |Me;p (see Fig. I14() . where blob represents the effective gluon 
propagator. 



B Collisional energy loss computations 

In this Appendix we will derive the collisional energy formula given by Eq. dJ. To do that we start 
from the Eq. (fT^ . i.e. 

spins 



where neq{k) is the equilibrium momentum distribution |24| at temperature T including quarks and 
gluons 

neq{k) = + J^^ (30) 

e|k|/T _ I g|k|/T _^ ]^ 

Here is the number of colors and Nf is the number of flavors. 

The matrix element M (see Eq. ((HI)), has been already computed in [25 (see Eqs. (45-46) in j24j ) . 
However, for completeness and due to a typographical error in [211, here repeat the main steps in 
the computation of \ Y^.pi^, \M\'^ . 



M = g^D^,{q)u{p', s')ru{p, s)u{k' , X')Yu{k, A). (31) 

In Coulomb gauge, the only non-zero terms in the effective gluon propagator are given in Eqs. @ 
and 0, which together with Eqs. (|ll3j) reduce the Eq. (PTT|) to 



M = g^/^L{q)u{p',s')fu{p,s)u{k',X')fu{k,\) 

+ g''^T{qW^ -ctq')u{p\s')^'u{p,s)u{k\}^)^^u{kA)- (32) 

Here = gVI^I- 

The matrix element given in Eq. (|32|) has to be squared, averaged over initial spin s of the jet and 
summed over all other spins. After evaluating the Dirac traces, and applying the assumption that 
|q| <C -E (highly energetic jet) we obtain similarly to 21| 
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spins 



+ 2Re{ALiq)AT{qy 



Iklfv . k' - "ll'^S + |k'|(v . k - ^4^) 



+ |Ar(g)P 



2 (v-k 



V • q q • k -.^ V • q q • k' 



|q| 



)(v-k' 



+ (|k||k'| -k-k')(u 



2 v-qq-v 



(33) 



In a static medium, the coUisional energy loss does not depend on the direction of v. Therefore, the 
Eq. H29|) can be further simphfied by averaging the integrand over the directions of v. The integrals 
that are required are 



Ji 



dn sin[(a;- v-q)4]^ 
47r (uj — V ■ q)2 



L 



4|q|i;2 
1 

4T;|q| 



Si{{v\q\ +uj)-)+ Si{{v\q\ - oj)-) 

V V 



1 - cos((i;|q| - a;)§) 1 - cos((u|q| + w)^) 



f |q| — (jj 



+ 



v\q\ + uj 



(34) 



and 



J2 = 



dn sin [(a;- v-q)|;]^ 
(w — V • q)2 



4tt 



(u - V • q) 



4f |q| 



V q - w - - C^ (u q + w -) + In 

f v|q| — UJ 



df] sin[(a;-v-q)^]2 ^ 2 1 A 
2 (t^ - V • q) = 1 



cos( — ) sm 



47r (u; — V • q^ 
By using Eqs. ()34l36p . it can be shown that 



(^|q|) \ 
^|q| ) 



(35) 



(36) 



dVL sin[(a;-v-q)|;]^ . . uj ^ q 

47r (w — V • qj^ |q 



= JiSt—r - JsT^ (37) 



and 

.7.^^iLl2 /„2|/J|2 2 



By using Eqs. (|33I38|) it is straightforward to show that 

d9. sin[(a;- v-q)^l2 



J 47r (lj — V • qj^ 



16,'i;'|A.to)|' '^l''l+f -I'^'V ., (39) 
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dn sin[(^-v-q)^ ]2 ^4 ^2 
An (a; — V • q)^ 



32/ i?2^e(AL(g)Ar((7)*) 



1,-,.^ v qq k' - v qq k 
|k|(v • k' ) + |k'|(v • k ) 



|q| 



(40) 



and 



dn sin[(a;- v-q)^]^ 
47r (cj — V • q)2 



16/ E'\ATiq)\' 



2 (V • k - ^4|^)(v • k' - + (Ikllk'l - k . ^')iv^ - 



16/i?^|Ar(g)| 



qr |q 

2(|q|2-^2)((2|k|+a;)2 + |q|2) ,,^,2 ,,2 



|qP 



|At(<?)P 



4|q|' 



[(t;2|q|2-a;2)Ji+2u;J2- J3] • (41) 



Therefore, averaging the f l-^P over the directions of v lead to 



'sin[(^-v.q)|;]2 1 ^ 



2\ 4z.2„A ,„„2(2|k|+^)2-|q|2 



16g''E'{\AL{qW 



(w - V • q)2 2 

spins 

2(|q|2-a;2)((2|k|+a;)2 + |q|2) ^ ,2 



Ji 



+|At(<z)|^ 



4|q|^ 



[{v^\q\^-oj^)Ji + 2ujJ2-J3]) 



(42) 



Since the colHsional energy loss does not depend on the direction of v, the Eq. ()29() can be written 



as 



AE, 



el 



Cr 



1 



^2 J (27r)32A; 
Cr 1 



neq[k) I ,^ 



(27r)32A;' 



'sin[(a;-v-q)|;]2 1 ^ 



(w - V • q)2 2 



spins 



327r4 ^2 



k||kV|k|rf|kVcosgne,(|k|)( r?^/"^' ^ E l-^l') ' (43) 



(w - V • q)2 2 



spins 



where 9 is the angle between vectors k and k'. Using the fact that q = k' — k, we obtain that the cos 9 
satisfies the following relation 



cos 9 = 1 



q2 — Ci;2 
2|k||k'| ' 



(44) 



where |k'| = |k| +u>. We can now introduce the change of variables from |k|, |k'| and cos0, to |k|, lo 
and |q|, which reduces the Eq. H29() to the following form 



AE, 



el 



Cr 1 
327r4 ^2 

CRg' 



^eq 

(|k|)(i|k| |q|ti|q| ojduj 



'sin[(w- v-q)^]2 1 



(cj - V • q)2 2 

spins 
\2 1^12 



neg(|k|)d|k| |q|d|q| tudu; (|AL(g)| 



2(2|k|+o;)2-|q|^ 



|2 , ,2\/^n|C| I , ,\2 I |r.|2N 



■Jl 



+ |Ar(g)| 



2(|q|2-^2)((2|k|+^)2 + |q|2) 



4|q|^ 



[(?;2|q|2 - c^2) + 2a; ^2 - J3] ) ■ (45) 
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Limits of the integration can be obtained from Eq. (|44() . from which it follows that 

< f Z < 2 , (46) 

2|k|(|k|+a;) 

leading to the limits in energy transfer lo 

Max[-|q|, |q| - 2|k|] < a; < |q|. (47) 

The limits on the momentum transfer |q| from coUisional scattering off a thermal parton with 
energy |k| is (see fTjl ) 

0<|q|<Min[i.,«^€^]. (48) 
l-t; + 2|k|/^ 

Here E and v are the energy and velocity of the jet. 

By using relations and (|15|) . the Eq. finally reduces to 

^Eel = neg(|k|)d|k| [j^ |q|d|q| j ^ Lodio + |q|d|q| j udco ] 



|q| J|k| ^|q|-2|k| 

2|k| 



(49) 

where |q| max is given in Eq. (|48|) . 
B.l Large L limit 

In this subsection we will consider the large L limit case, and compute the coUisional energy loss per 
unit length. To do that we multiply both sides of the Eq. by ^ i.e. 

— A^eZ = J n,,{\k\)d\k\ |q|d|q| udu (|Ai(g)P^ ^ ^^^1 

W^^^K^^ - .^)^J. + 2.^ J. - ^^3)). (50) 

4|qp ttL ttL ttL 



To compute :^c7i,2,3 in the limit when L — > oo we will use the following expression 

2z; sin[(a;-v-q)|;]2 ^ 

- >o(t(; — v-q). (51) 



vrL (u; — V • q) 



Then, 



2v L— oo 




/^^(.-v.q) = ^e(.W-^), (52) 



2?; L^oo f dQ 
—rJ2 — 



y"^5(cu-v-q)(a;-v-q) =0 (53) 



and 
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^^3^ I ^5(t.-v.q)(u;-v.q)2=0, (54) 

leading to 

^AE,i = £g!yneg(|k|)d|k| |q|d|q|u;da;^e(7;2|q|2-^2) 

Therefore, by using relations 1)47^ . ()48() and i^'^jqp > w^, the coUisional energy loss per unit length 
in an infinite size QCD medium reduces the following expression (Cr = 4/3) 

—TT = a 2 3 / neq{\\<i\)d\\<i\ / d|q| / ujduj + d\q\ / ^ wdo; 

dL 6v^TT-^Jq J|k|/(l+i,) -'|q|-2|k| / 

/ ,„ (2|k|+u;)^-|qp , ,^ ^„ (|qp-a;^)((2|k|+^)^ + |qp) ^ \ 
M^lWI ^ + |Ar(g)| ^-^4 [v |q| -w ) I . (56) 
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